Galaxy rotation curves from a fourth order gravity 



Priti Mishra and Tejinder P Singh 

Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai, India 
E-mail: priti@tif r .res . in, tpsingh@tif r . res . in 

Abstract. While the standard and most popular explanation for the flatness of galaxy rotation 
curves is dark matter, one cannot at this stage rule out an explanation based on a modified law 
of gravitation, which agrees with Newtonian gravitation on the scale of the solar system, but 
differs from it on larger length scales. Examples include Modfied Newtonian Dynamics [MOND] 
and Scalar- Tensor- Vector Gravity [STVG]. Here we report on a fourth order modification of the 
Poisson equation which yields the same Yukawa type modification of Newtonian gravity as 
STVG, and which can explain flat galaxy rotation curves for a large sample of galaxies, once 
specific values for two parameters have been chosen. We speculate on two possible origins for 
this modified Poisson equation: first, a possible fourth order modification of general relativity, 
and second, quadrupole gravitational polarization induced on a galaxy because of the pull of 
neighbouring galaxies. 



1. Introduction 

Flat galaxy rotation curves and the departure from the expected Keplerian fall-off which they 
imply are most popularly attributed to the presence of cold dark matter [CDM] in the galaxy. 
Together with a cosmological constant A, CDM constitutes the highly successful ACDM model 
which explains most of cosmological, cluster and galaxy data. Nonetheless, considering that 
dark matter has not yet been detected in the laboratory, one should tentatively also allow for 
the possibility that flat galaxy rotation curves arise, not from the presence of dark matter, but 
from a possible modification of the law of gravitation on scales larger than the solar system. 

Notable examples of such modifications include Modified Newtonian Dynamics [MOND] [1 
where it is assumed that the gravitational force acting on a test particle of mass m is given by 
the relation 

F = mau(—] (1) 
\ao/ 

where a is acceleration in Newtonian mechanics and ao is a new fundamental constant of nature 
having the value 2 x 10 -10 ms~ 2 . For very small accelerations [large distances] it is assumed 
that 
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whereas \x approaches the value one for accelerations encountered in the solar neighborhood. 
Hence for large distances we will have the relation 
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The virtue thus is that a falls as 1/r rather than 1/r 2 . Thus equating a to the centripetal 
acceleration v 2 jr we get that 
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With the numerical choice of ao made above one gets the desired flat value of the velocity. 

Another noted example is Scalar- Tensor- Vector Gravity [STVG], where a massive vector field 
with a mass, denoted fj,, interacts with gravity with a coupling strength denoted u> [2]. It can 
be shown that in STVG, the law of gravitation acquires a Yukawa modification, which is given 
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Here, r$ = 1/n and Mq is a parameter which vanishes when cj = 0. It is assumed [and justified 
by further considerations] that one can generalize this to the case of a mass distribution by 
replacing the factor GM/r 2 in ([5]) by GM(r)/r 2 . The rotational velocity of a star v c is obtained 
from v 2 (r)/r = a(r) and is given by 
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From observations of galaxies the density profile fj,(r) is known to be 
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1 for HSB galaxies, 

2 for LSB & Dwarf galaxies. 



A good fit to a large number of galaxies has been achieved with the parameters: 
M 



9.60 x 10 11 M & , r = 13.92 kpc = 4.30 x 10 22 cm. 
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In the fitting of the galaxy rotation curves for both LSB and HSB galaxies, using photometric 
data to determine the mass distribution M(r), only the mass-to- light ratio (M/L) is employed, 
once the values of Mq and ro are fixed universally for all LSB and HSB galaxies. Dwarf galaxies 
are also fitted with the parameters 
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2.40 x 10 11 Mq, r = 6.96 kpc = 2.15 x 10 22 cm. 



(11) 



2. Fourth order gravity 

Remarkably enough, it can be shown that the modified law of acceleration (0) of STVG, and the 
consequent velocity profile (0) is also a solution of the following modified fourth order Poisson 
equation [3] 

V 4 - k 2 V 2 4> = -4irGk 2 fi(r) (12) 

provided the density profile is chosen as in ([7]) and one identifies k = 1/vq. Mq appears as a 
constant of integration in the solution of (|12p . 



Greater insight into the modified acceleration can be had by expanding the acceleration law 

(13) 
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around r = tq and writing it as a sum of two terms : a part that falls as 1/r 2 and is independent 
of ro, and a part that falls as 1/r and depends on tq. This gives 
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Figure 1. Dashed curve: velocity due to first term, dotted curve: velocity due to second term, 
thick curve: total velocity. 



The second term dominates for r > r$. Fig. 1 above plots the rotation velocity curve, and 
it is clear that the Keplerian fall-off due to conventional Newtonian acceleration [first term] is 
modified to a flat rotation curve [second term] for r > ro. The 1/r fall-off brought about by the 
second term is MOND-like behaviour. At even larger distances, the second term is exponentially 
damped and the fall off is again Keplerian but with an effectively larger value of G. 

It is significant that from our work we can give an estimate of the theoretical value of ao in 
MOND. A simplistic guess would be to construct a quantity with the dimension of acceleration 
from our fundamental quantities k and Mq. This is nothing but 

6.67 x KT 11 x 9.60 x 10 11 x 2 x 10 30 

(13.92 x 10 3 x 3.08 x 10 16 ) 2 ^ ' 

3 x 10~ 10 ms~ 2 (16) 

which is of the same order as ao in MOND. ( ao = 2 x 10 _10 m/s 2 ). The relation is made more 
transparent by comparing the acceleration given by the second term of Eqn. (|14|) with the 
acceleration in MOND in Eqn. Q. This comparison yields the fundamental relation 
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3. Possible origin of fourth order gravity 

The fourth order modified Poisson equation considered above is the weak field limit of the 
following modification of general relativity: 
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It is possible that such an effective modification of Einstein gravity comes into play during the 
late stages of the evolution of the Universe, when large scale structures form. 

An intriguing alternate possibility, which has been discussed in [3] is that the fourth order 
modification of the Poisson equation is a consequence of induced polarization of the gravitational 
field of a galaxy, caused by neighboring galaxies. Thinking of galaxies as 'molecules' made of 
atoms [the stars], one would like to analyze how the averaged gravitational field inside a galaxy, 
is modified by the polarization of the molecules, due to the external pull of other galaxies. Indeed 
one has at the back of the mind the polarization induced modification of electromagnetic fields 
in an electrically charged medium. 

It can be shown that as a consequence of polarization, the averaged Einstein equations are 
effectively modified to 




where T^i is the energy-momentum tensor of molecules, and Q^ pU a is the quadrupole 
gravitational polarization tensor. It is modelled by assuming it be proportional to the Weyl 
tensor, in analogy with the electrostatic case, 



QiOjO — ZgCiOjO (20) 

and the gravitational dielectric constant e g is modelled by [I] as 
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Here, A is the typical linear dimension of a molecule, m is the typical mass of a molecule, ooq is a 
typical frequency of harmonically oscillating atoms in the molecule and N is the number density 
of the molecules in the medium. In the weak field limit, Eqn. (|19p can be shown to go over 
to the modified Poisson equation (|12p with k 2 = c 2 /2irGe g [5]. Furthermore, the theoretically 
predicted values of k and Mq agree well with the values phenomenologically assumed for fitting 
data, suggesting that gravitational quadrupole polarization is a contender for explaining flat 
galaxy rotation curves [3]. 
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